Splitting Sensitivity of the Ground and 7.6 eV Isomeric States of ^^^Th 
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The lowest-known excited state in nuclei is the 7.6 eV isomer of Th. This energy is within 
the range of laser-based investigations that could allow accurate measurements of possible temporal 
variation of this energy splitting. This in turn could probe temporal variation of the fine-structure 
constant or other parameters in the nuclear Hamiltonian. We investigate the sensitivity of this 
transition energy to these quantities. We find that the two states are predicted to have identical 
deformations and thus the same Coulomb energies within the accuracy of the model (viz., within 
roughly 30 keV). We therefore find no enhanced sensitivity to variation of the fine-structure constant. 
In the case of the strong interaction the energy splitting is found to have a complicated dependence 
on several parameters of the interaction, which makes an accurate prediction of sensitivity to tem- 
poral changes of fundamental constants problematical. Neither the strong- nor Coulomb-interaction 
contributions to the energy splitting of this doublet can be constrained within an accuracy better 
than a few tens of keV, so that only upper limits can be set on the possible sensitivity to temporal 
variations of the fundamental constants. 

PACS numbers: 23.20.-g,06.20. Jr,27.90.+bb,42.62.Fi,21.10.Sf,21.60.-n,21.60.Ev 



The isotope ^^^Th has recently become of interest be- 
cause of its unusually low- lying (7.6 eV, 3/2+) [l| first ex- 
cited state, which is an isomeric state with an estimated 
half-life of 5 hours. The 5/2+ — 3/2+ ground-state-to- 
isomer transition energy is within the range of atomic 
transitions, and it has been suggested that this almost 
degenerate doublet in ^^^Th could be used as a nuclear 
clock Q and as a sensitive probe of possible temporal 
variation of fundamental constants, including the fine- 
structure constant (a) and the quark mass^Sj. The sen- 
sitivity of the transition energy to temporal changes in 
the fundamental constants varies considerably depend- 
ing on the assumptions madefsl, 0, Q. For example, in 
refs. 0, Q temporal variation of the fine structure con- 
stant (d) was related to constants in the Nilsson Hamil- 
tonian, whereas in [1] d was shown to be proportional to 
the Coulomb energy difference between the two states, 
which requires detailed information about the deforma- 
tions involved. Thus it is important to understand the 
nuclear-structure issues giving rise to this doublet. 

In the present work we examined this doublet using the 
finite-range microscopic-macroscopic model (FRDM) 6], 
which describes many nuclear-structure properties (such 
as ground-state masses and deformations) over a broad 
range of nuclei. Our goal is to examine the sensitivity 
of the energy splitting between these states to the un- 
derlying components of the effective nuclear interaction, 
including the single-particle potential and the pairing, 
spin-orbit, and Coulomb interactions. Knowledge of this 
sensitivity is essential for determining the sensitivity of 
the transition energy to possible temporal variation of 
fundamental constants. 

In our macroscopic-microscopic model(6i] the macro- 
scopic terms give the smooth variation of the nuclear po- 
tential energy (mass) with proton number Z, neutron 
number N , and deformation. The dependence of nuclear 
structure properties on microscopic quantum-mechanical 
effects is obtained from a deformed single-particle poten- 



tial through the use of Strutinsky's method0, Nine 
constants of the macroscopic part have been determined 
in a least-squares adjustment to 1654 measured nuclear 
masses with Z > 8 and > 8; the details are given in 
0. For these 1654 nuclei ranging from ^^O to ^^^^106 
the model mass accuracy is 0.669 MeV. For the heavier 
mass regions the model is more accurate: for N>65 the 
corresponding accuracy is 0.448 MeV. Values of other 
model constants (such as the depth and diffuseness of the 
single-particle potential, and the strengths of the spin- 
orbit and pairing interactions) are determined from other 
global considerations, as discussed in Ground-state 
masses and shapes have been calculated for 8979 nuclei 
and tabulated in [6] . The shape parameters tabulated are 
quadrupole (£2), octupole (ea), hexadecapole (£4), and 
hexacontatetrapole (ee) deformation (shape) degrees of 
freedom. Strong evidence of reliability of the model is 
given by its now well-established predictive capabilities 
for new nuclear-mass regions, including unstable nuclei 
and super-heavy nuclei. 

In the present calculations we carried out a high- 
accuracy determination of the ground-state deformation 
and quasi-particle energy by minimizing the potential en- 
ergy for ^^^Th on a fine deformation grid. Taking ad- 
vantage of enhanced present-day computational power, 
we varied all four shape parameters in steps of 0.001, 
which is to be compared to the 0.05 grid step used in 0]. 
We predict three very closely lying states with identical 
deformation: the 5/2+ ground state, the 3/2+ isomer, 
and a 5/2^ state. The three lowest (almost degenerate) 
neutron quasi-particle states are predicted to have de- 
formations: £2 = 0.170, eg = 0.000, 64 = -0.084, and 
eg = —0.002. The 5/2^ state, with asymptotic quan- 
tum numbers [iVn^A, if'^] — [732,5/2^], lies lowest in 
energy. The ground-state and isomer doublet of inter- 
est have asymptotic quantum numbers [6 3 3,5/2+] and 
[6 31, 3/2+], respectively, and are predicted to he at 13.6 
keV and 21.6 keV above the 5/2". 
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Since the three lowest neutron quasi-particle configura- 
tions are predicted to have the same deformation (which 
determines their proton distribution) , they also have the 
same Coulomb energy (within the model uncertainties). 
This is not the case for other states in ^^^Th, for which 
the energy was found to be minimized by different values 
of the shape parameters. 

Although minor changes to the model parameters 
could move this triplet of states relative to one another, 
the exact position of the 5/2~ state is not essential to 
our present interest, since the predicted energy of that 
state does not directly affect the sensitivity of the split- 
ting between the 5/2+ — 3/2+ doublet to the underlying 
interaction. Given the global accuracy of the model in 
this mass region, we expect that the 5/2~ state should 
lie somewhere within 50 keV of the ground state. 

The calculated positions of the three lowest-lying 
states as functions of deformation are shown in Fig. 1, 
and the three states are seen to track very closely in en- 
ergy for values of €2 between roughly 0.12 and 0.22. As 
discussed below, we find a similar close tracking for all of 
the deformation parameters: £2, • • • , ee- For comparison, 
we also show in Fig. 2 the [6 31,3/2+] and [6 3 3,5/2+] 
doublet together with two other low-lying states of ^^^Th. 
Note the very different energy scales used in Figs. 1 and 
2, and that the latter two states do not show a minimum 
at the same deformation as the three lowest states, which 
is the more typical situation seen in deformed nuclei. The 
energies of different quasi-particlc configurations for a nu- 
cleus are not generally minimized by the same values of 
deformation parameters. 

The predicted energy separation of the [6 31, 3/2+] and 
[6 3 3,5/2+] doublet (viz., 8.3 keV) is large compared to 
the actual excitation energy of the isomer (7.6 eV), but 
it nonetheless means that on a "normal" nuclear-energy 
scale the two states are predicted to be almost degener- 
ate. For comparison, typical quasi-single-particle energy 
splittings in this mass region are a few hundred keV. For 
completeness we note that the odd-neutron wave func- 
tion for the ^l^ = 3/2+ isomer is calculated to have the 
following asymptotic Nilsson components: 

|3/2+) = 0.650[633] + 0.519[642] + 0.431[613] 

-G.221[622] - G.137[6G2] - 0.131[853] , (1) 

while that of the fl'^ — 5/2+ ground state has compo- 
nents: 

|5/2+) = 0.656[631] + 0.521[642] + 0.403[611] 

-0.242[651] - 0.123[871] - 0.104[622] . (2) 

Those components in the wave functions whose squared 
amplitudes are less than 1% are not listed. 

When quasi-particle configurations have the same de- 
formation, the energy difference between them is given 
by 

Ai?g.,._i3o = v/(eiso - A)2 + A2 - ^(eg.,. - A)2 + ,(3) 



where A is the Fermi energy, A is the pairing gap obtained 
by solving the pairing equations, and are the single- 
particle energies for the two states; all of these energies 
depend on deformation. The near degeneracy between 
the ground state and isomer of ^^^Th arises because the 
Fermi surface energy is about midway between ejso and 
Eg.s. . The isomer and the ground state will be exactly 
degenerate when (ciso — A) = — (cg.s. — A) (i.e., for A = 
(fiso + eg.s.)/2). 

The single-particle energies Si, depend on the deformed 
single-particle potential well and on the spin-orbit in- 
teraction. These in turn depend on parameters govern- 
ing the depth and diffuseness of the potential well and 
the strength parameter of the spin-orbit interaction. Al- 
though the single-particle energies depend on the shape 
of the potential, the shape parameters are not in the cat- 
egory of adjustable parameters, and single-particle ener- 
gies are calculated by minimizing the energy of the chosen 
configuration. 

The dependence of the doublet splitting on the fine- 
structure constant is particularly straightforward to ex- 
amine in our model, as we discuss below. On the other 
hand, expressing that splitting in terms of more fun- 
damental parameters of the strong interaction (such as 
the quark mass, for example) would require a very de- 
tailed and non-trivial analysis of the relation between the 
model parameters and fundamental sub-nucleon degrees 
of freedomQ. Our analysis is therefore necessarily re- 
stricted to a study of the sensitivity of the predicted dou- 
blet splitting to the effective interactions in our model. 

The similarity in the shape of the ground and isomeric 
states of ^^^Th implies that (within the accuracy of the 
model) they have the same charge (i.e., proton) distri- 
bution and therefore the same Coulomb energy. This 
has unfortunate implications with respect to experimen- 
tal searches for a temporal variation in the fine-structure 
constant (d) obtained from measuring a temporal vari- 
ation in the energy splitting between these two states 
(denoted by lu). The latter variation is proportional to 
the Coulomb-energy difference of the two states [4] : 

^=((^c))-, (4) 
a 

where ((Vc)) = (iso | Vc | iso) - (g.s. | Vc | g.s.). 
We therefore find essentially no sensitivity to a, since 
((Vc)) = for our chosen mesh-parameter step size. Al- 
lowing £2 to vary by the mesh-parameter step size (0.001) 
in this calculation produces a variation in ((Vc)) of ap- 
proximately 30 keV, and this sets the upper limit for 
variations in to relative to a/ a. Uncertainties in the ef- 
fective nuclear interaction prevent any nuclear-structure 
calculation from being predictive on an eV scale. 

With the exception of accidental degeneracies, nearly 
degenerate doublets in deformed nuclei reflect a strong 
similarity in the deformation of the states involved. If 
the deformations of two states are not similar their en- 
ergy splittings are typically at least several tens of keV. 
To illustrate this and to investigate the accuracy of the 
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^^/g.S. — ISO 


Best 


v;,(+io%,-io%) 


a(+10%,-10%) 


A„(+10%,-10%) 


G(+10%,-10%) 


8.8 keV 


-4.4keV, 31.6 keV 


25.4 keV, 2.0 keV 


-44.2 keV, 29 keV 


8.7 keV, 9.2 keV 



TABLE I: Variation in the predicted 3/2"'"-5/2"'" doublet energy splitting for induced by a ±10% variation of four of the 

global parameters of the model (see Ref.[^). The column labeled "Best" is the predicted splitting with the standard values of 
these parameters, Vb is the depth of the single-particle potential, a is the range or diffuseness of the single-particle potential, 
A„ is the strength of the neutron spin-orbit interaction, and G is the pairing strength. Splittings that are listed with a negative 
sign mean that the [6 3 3,5/2^] and [6 31,3/2^] states were inverted. 
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FIG. 1: [Color online] Calculated quasi-particle energies 
for the three lowest quasiparticle states (relative to the 
[6 3 3,5/2+] ground state, which therefore has zero energy for 
all deformations) in ^^^Th as functions of e2. All three states 
remain almost degenerate at all values of £2 near the common 
minimum at £2 = 0.170. 



FIG. 2: [Color online] Calculated quasi-particle energies of 
the 3/2+-5/2+ doublet and of a 7/2" state and a 13/2+ state 
in ^^^Th as functions of £2. The doublet states remain almost 
degenerate for a wide range of values of £2. Note the change 
in energy scale from Fig 1. 



model for very close-lying (eV) doublets we examined the 
ground and isomeric states in ^^^U, where the observed 
energy splitting is 76 eV. We calculated this (7/2~, 1/2+) 
doublet, again allowing all of the deformation parame- 
ters to vary independently for each state. The results are 
very similar to those for ^^^Th, with the ground-state 
and isomer energies being minimized by the same values 
of the deformation parameters, and lying just below and 
above the Fermi surface, respectively. The deformations 



for these two states in ^ssy -v^ere found to be £2 = 0.205, 
£3 = 0, £4 = —0.07 and £5 = 0.025, and the predicted 
energy splitting was 29 keV. 

We next considered the sensitivity of the structure of 
the ^^^Th doublet to several of the globally determined 
strong-interaction parameters of the model, namely, the 
depth and diffuseness of the single-particle potential, 
and the strength of the pairing and spin-orbit interac- 
tions. We varied each of these four global parameters by 
±10%. We note, however, that variations of this mag- 
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nitude would likely destroy the agreement with ground- 
state masses and for some nuclei would lead to incorrect 
assignments. For example, the strengths of the neutron 
and proton spin-orbit interactions have been fitted to ex- 
perimental levels in the rare-earth and actinide regions, 
where properties of these nuclei (particularly level spac- 
ings and orderings) strongly suggest a linear dependence 
of the spin-orbit strength on the nuclear mass number A: 

An,p — ^n,p A -\- ln,p • (5) 

With this parameterization the constraints on 
fcn, fcp, ^n, ^p are considerably tighter than 10%. Nonethe- 
less, in order to gain a better understanding of the origin 
of the near degeneracy in ^^^Th we varied the strength 
of this and the other three parameters by ±10%. The 
results are summarized in Table 1. We find that the 
predicted doublet splitting varies on the keV scale 
within a 10% variation of these four global parameters, 
while for some variations the levels arc predicted to be 
inverted (indicted by a mirnis sign in Table I). There is 
no one parameter or part of the interaction determining 
or dominating the magnitude of the predicted splitting. 
The degeneracy instead arises from a complicated 
combination of different effects causing the ground state 



and isomer to lie just below and above the Fermi surface, 
respectively. 

In summary, we have examined the structure of the 

ground state and isomer in ^^^Th. Within the present 
model the two states are predicted to be different quasi- 
single-particle neutron states, but corresponding to the 
same nuclear deformation. The Coulomb energy is pre- 
dicted to be the same for the two states (within an uncer- 
tainty of roughly 30 keV induced by our mesh-step size) , 
which suggests no enhanced sensitivity to a. The origin 
of the near degeneracy appears to be accidental and is 
difficult to parameterize in terms of any one component 
of the effective interaction. Nevertheless, this doublet re- 
mains of interest in possible searches for time variations 
in fundamental physics because of the small energy split- 
ting involved. 
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